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ABSTRACT 

In the widely adopted ACDM scenario for galaxy formation, dwarf galaxies are the building blocks 
of larger galaxies. Since they formed at relatively early epochs when the background density was 
relatively high, they are expected to retain their integrity as satellite galaxies when they merge to 
form larger entities. Although many dwarf spheroidal galaxies (dSphs) are found in the galactic halo 
around the Milky Way, their phase space density (or velocity dispersion) appears to be significantly 
smaller than that expected for satellite dwarf galaxies in the ACDM scenario. In order to account 
for this discrepancy, we consider the possibility that they may have lost a significant fraction of their 
baryonic matter content during the first infall at the Hubble expansion turnaround. Such mass loss 
arises naturally due to the feedback by relatively massive stars which formed in their centers briefly 
before the maximum contraction. Through a series of N-body simulations, we show that the timely 
loss of a significant fraction of the dSphs initial baryonic matter content can have profound effects 
on their asymptotic half-mass radius, velocity dispersion, phase-space density, and the mass fraction 
between residual baryonic and dark matter. 

Subject headings: Methods: numerical, Galaxies: dwarf. Galaxies: evolution. Galaxies: structure, 
(Cosmology: ) dark matter 



1. INTRODUCTION 

Dwarf galaxies play an important role in the ACDM 
scenario for galaxy formation. They form early, in rela- 
tively dense background environments, and are the build- 
ing blocks of larger galaxies. ACDM simulations show 
that these dwarf galaxies retain their integrity during 
their merger process and predict a rich population of 
satellite dwarf gala xies around large galaxies such as 
the M ilky Way (e.g. iDiemand et"am2008l : ISpringel et all 

Hool). 

Many satellite dwarf spheroidal galaxies (dSphs) in 
the Galactic halo have been found through SDSS and 
other surveys. Follow-up spectroscopic observations re- 
veal these dSphs are surrounded by dark-matter ha- 
los. But, the internal phase space density of the 
most newly discovered (faint) satellite dSphs is sig- 
nificantly smaller than that predicted from the dark- 
matter-only simulations. This discrepancy poses a 
challenge to t he standard ACDM scenario in its sim- 
plest f orm (e.g.lBovlan-Kolchin et al.l[2"011l:lFerrero et all 



20111: iRashko v et al.l 120121: iBovlan-Kolchin et all 120121: 
Wolf fc Bullocy|20li). Recent observations show that 
this problem is not limited to th e Milky Way, bu t present 
in Andromeda as well (e.g. .ToUerud et al.ll2012l ). 

We examine the dSphs' structural adjustment associ- 
ated with the loss of baryonic matter due to star for- 
mation feedback. Today, the dSphs' total dark to bary- 
onic matter ratio is generally much larger than the cor- 
responding cosmological value (jMatco 1998). The pref- 
erential loss of baryonic matter may be the potential so- 
lution to both the "mi ssing baryonic matter" and the 
missmg satellite" (e.g iKlvpin et al.l 119991: iMoore et all 
Il999f) puzzles. 
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There is very little detectable molecular or atomic gas 
inside the dSphs today. In ^ we briefly discuss the 
possibility that a significant fraction of the initial gas 
content may be lost due to the feedback effect of first 
generation stars in them. Presumable, this effect is most 
intense during the collapse following the dSph progeni- 
tors' turn around from the Hubble flow. We show, in ^J3l 
how the loss of a small fraction of the total mass may sig- 
nificantly affect dSphs' asymptotic dynamical structure. 
Baryonic matter may also be removed from the dSphs 
through tidal or r am pressure strippin g by tenuous gas in 
the Galactic halo (|Lin fc Fabeijll983| ). However, gradual 
loss of gas from the dSphs' virialized potential would not 
significantly affect their internal dynamical structure. 

There have been several attempts to include bary- 
onic matter in comp rehensive cosmologi cal simulations 
of galaxy formation ()Gnedin et al.l [20091 ). For compu- 
tational efficiency, a simple-to-use prescription for star 
formation rates (as a function of local gravity, gas den- 
sity, temperature, metallicity, and background radiation) 
is clearly technically desirable. However, there are many 
competing effects including atomic processes, radiation 
transfer, hydrodynamics and magnetic fields, which are 
sensitive to the diverse initial and boundary conditions 
over a large dynamic range in spacial extent and time 
scales. The determining factors for the onset, efficiency, 
and impact of star formation remain enigmatic. 

If the present-day stellar content is distributed widely 
in the form of initial gas through dSphs, the background 
UV flux after the epoch of re-ionization would easily 
ionize the gas within typical dwarf galaxies, delaying 
or ev en preventing cooling and star form ation within 
them (|Dong et al.ll2003l : iBrown et al.ll2012D . Most dSphs 
within the Local Group are observed to contain multiple 
generations of stars, the oldest of which formed in the 



2 



Matthias Gritschneder, Douglas N.C. Lin 



early epochs of cosmic evolution, when the background 
UV flux may have been intense. Possible solutions to this 
conundrum include that the formation of these stars pre- 
ceded the epoch of re-ionization or that these dSphs may 
have lost most (up to ^ 90%) of their original baryonic 
content which provided a more effective shield against ex- 
ternal UV photons. The second possibility is in line with 
the large mass ratio between dark and baryonic matter 
in dSphs. 

Most recent cosmological simulations 

(jBrooks fc Zolotovll2012tlGovernato et al.ll2012.1 produce 
results which appear to be in agreement with our 
baryonic mass loss conjecture. It has also been argued 
that about 10% of the haloes in sim ulations might not 
face the 'too big too fail' problem (jPurcell fc Zentneil 
I2OI2I ). However, this conjecture only shifts the problem 
at hand towards the uniqueness of the Milky Way. Here, 
we present an idealized toy model for the formation of 
dSphs. We introduce a prescription to link the loss of 
the interstellar medium with the burst of first generation 
star formation in dSphs. The main advantage of our 
approach here is that it enables us to identify the key 
physical effects which determine dSphs' asymptotic 
dynamical structure. On the technical side, it also 
provides adequate resolution on the scale of the dSphs. 

We first review, in ^ the basic assumptions, method- 
ology, initial and boundary conditions used in our nu- 
merical models. We then continue to present the results 
of our dimensionless simulations (*|3|). In §4] we convert 
these numbers to physical quantities and in f}5]we draw 
the conclusions. 

2. BASIC APPROACH AND INITIAL CONDITIONS 

The main cause of the loss of ordinary matter may 
be feedback associated with s tar formation. In previous 
simulations (jPong et a l.||2003f) on the effects of radiative 
transfer and photoionization in dSphs, it was shown that 
star formation can indeed proceed in their more massive 
and gas rich progenitors. The star formation and gas re- 
tention efficiency may vary widely in galaxies with simi- 
lar dark matter potentials because they depend on many 
factors, such as the baryonic fraction, external perturba- 
tion, initial mass function, and background UV intensity. 
The presence of very old stars in dSphs indicates that 
their initial baryonic-to-dark matter content was com- 
parable to the cosmic value. This constraint suggests 
that the initial density fiuctuation of baryonic matter 
may be correlated with that of the dark matter. For the 
more massive dwarf elliptical galaxies, the star forma- 
tion efficiency and gas retention rate are much higher. 
Their mass-to-light ratio is regulated by star formation 
feedback and is expected to be nearly independent of 
their absolute luminosity. The results of our theoretical 
models repro duce the observed {M/L) — correlation 
(|Mateolll998h . 

Many dwarf galaxies contain stars with diverse [Fe/H] 
and some contain multiple generations of stars, includ- 

ing a recent Gyr) episode in the Fornax dSph 

(|de Boer et all l2012[ ). Significant iron dispersion in 
dSphs suggests they were self contaminated by the re- 
tained and recycled supernova ejecta from massive stars. 
But, the total metal content in these systems is much 
less than expected from the heavy-element production of 
massive stars in each episode of star formation. Such 



a deficiency implies that a substantial fraction of the 
dSphs' gravitational potential is generally shallow with 
a modest velocity dispersion. The UV flux from a few 
massive supernova- progenitor stars would be adequate 
to ionize the residual ga s and evaporate it from the dSphs 
(jNoriega-Crespo et al.llT989i ). In the presence of an ini- 
tially disturbed gas distribution (due to photoionization) , 
only a few percent of the gas en riched by the supernovae 
remains m the center of dSphs (jFragile et al.l [2003'). 

The ionization and supernova blast wave only directly 
affect the dynamics of baryonic matter. Although the 
ejection speed may substantially exceed the escape speed 
of the dSphs' dark matter halo, the loss of baryonic mat- 
ter makes small fractional changes in their total mass and 
potential. If this mass loss occurs after the dark matter 
halo has already virialized, it will only modify the dark 
matter content by a small amount. However, we are able 
to show here that a timely loss of baryonic matter during 
the virialization process would lead to an extended halo 
structure which is prone to subsequent harassment. 

During the initial infall of the dSphs, both dark and 
baryonic matter attain infall velocities comparable to 
their escape speed before virialization. In the absence 
of any feedback, the gravity of the dark matter potential 
accelerates both components to attain a kinetic energy 
which is slightly less than their gravitational energy. As 
the density of both dark and baryonic matter increases, 
gas in the dSphs become self shielded against the back- 
ground UV radiation. Thermal instability leads to the 
formation of dense cores and young stars. Photoioniza- 
tion and supernovae feedback induces the ejection of gas 
from the central regions of dSphs. At this advanced stage 
of infall, the total energy of the infalling dark and bary- 
onic matter is a small fraction of the instantaneous po- 
tential energy (lAarseth et al.|[l988t ) and a small fractional 
mass loss could offset the near cancellation between the 
kinetic and potential energy. 

We consider the effect of baryonic matter loss on the 
evolution of the dark matter halo in dSphs, especially 
during the infall stage of their formation. In order 
to illustrate this effect, we consider a series of N-body 
simulations with th e publicly available SEREN code 
(jHubber et al.l[20Tll ). The dark and baryonic matter is 
represented by 10000 particles which only interact with 
each other through their mutual gravity. For computa- 
tional simplicity, we represent the Hubble expansion turn 
around with systems of particles at rest with a uniform 
initial density. 

The system contracts under its own gravity. By re- 
ducing the mass of individual particles, we impose an 
instantaneous loss of a prescribed fraction / of the total 
mass after the half mass radius i?hm of the systems has 
reduced to i?ioss from their original total radius Rq by a 
factor up to 20 (see Sj3|). Our goal is to determine the 
asymptotic structure of the systems after virialization as 
a function of / and -Rioss- Therefore, we need to generate 
models with the minimum value of i?hm < R\oas- 

The minimum value of i?hm/^o is determined by the 
initial fiuctuation spectrum. For systems represented by 
N particle s, its value, due to th e global random noise, is 
(jAarseth et al.lll988[ ). This contraction factor 
would be limited in cosmological simulations with inad- 
equate resolution on the scale of dSphs. For the entire 
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Fig. 1. — Dependence of the half-mass radius on the precise time 
of loss for different times of loss, characterized by the contraction 
of the halo at the time of the loss. We plot the radius in units of 
the initial radius and the time in units of the free-fall time. 
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Fig. 2. — Dependence of the half-mass radius on the amount of 
loss for cases with a loss between 5% to 25% of mass loss. We plot 
the radius in units of the initial radius and the time in units of the 
free-fall time. 

range of models in Jj3l it is adequate to use N = 10'*. Al- 
though the minimum value of Rhm/ Ro may reduce with 
larger values of N, the asymptotic structure is affected by 
the global changes in the gravitational potential rather 
than two body relaxation process (i.e. it is determined 
by the values of / and R/Rq rather than N). 

For a fiducial model, we neglect any mass loss. In this 
case, the system relaxes into a virial equilibrium after 
a few times of its initial free fall time scale tg. A small 
fractional mass loss, especially during the advanced stage 
of the infall, can lead to much more extended relaxed 
kinematic structure after violent relaxation. We run the 
model for 2te, which is sufRcient to reach a stable half 
mass radius in the fiducial case. 

3. RESULTS 

Here, we investigate a number of different scenarios. 
We vary the precise time of the mass expulsion, the 
amount of mass ejected as well as the region where the 
mass is expelled. 

First, we investigate the timing of the loss. We investi- 
gate cases with a contraction to 1/2, 1/4, 1/10 and 1/20 



Fig. 3. — Dependence of the half-mass radius on the region of 
loss. 10% of the material is lost in to within a range of different 
regions. We plot the radius in units of the initial radius and the 
time in units of the free-fall time. 

at the epoch of the mass loss, respectively. The cor- 
responding dynamical timescales (to bounce back) are 
O.aStff, O.iatff, O.Oatff, and O.OUff. From Figure ffl it be- 
comes clear the timing is of crucial importance. The 
closer the loss to the maximum contraction, the more 
drastic the effect of baryonic mass loss on the entire 
dark matter halo. From a physical point of view, the 
first burst of star formation and thus the most efficient 
feedback is indeed to be expected at the first time high 
densities are reached, i.e. close to maximum contraction. 

Another interesting feature is the amount of mass loss 
needed. Clearly, this parameter is strongly constrained 
by cosmological parameters, i.e. the maximum mass loss 
is limited by the amo unt of baryonic matter initially. Ac- 
cording to WMAP7 ()Komatsu et al.ll201lD . the ratio of 
baryonic to dark matter is about 1:4.45. Assuming this 
ratio is preserved in the initial halo, there is about 17% 
baryons inside the halo. We investigate a series of cases 
with an loss of 5% to 25%. Above 25% mass loss the 
halo disperses. Figure [2] shows a linear dependence of 
the final half-mass radius on the amou nt of mass loss (as 
expected from e.g. lAarseth et al.lll988l ). An encouraging 
result is that even the loss of about 10% of the total mass 
(i.e. slightly more than half of the baryonic mass) leads 
to a significant puff up, resulting in an halo twice a big 
as in the fiducial case. 

In the context of the observed (M/L) — Mv relation 
and the fi nal ratio of baryonic to d ark matter of about 
1:100 fe.g. lJardel fc GebhardtllMl) . an even higher frac- 
tional mass loss (~ 15%) is desirable. However, here we 
consider only a loss of 10% during the first infall, which 
is sufficient to explain the low phase-space density. The 
observed mass to light ratio is most likely later achieved 
via star bursts or tidal effects. 

To increase the applicability of this numerical experi- 
ment to the problem at hand, it is important to constrain 
the region where the mass loss occurs. In a dwarf galaxy, 
the mass loss is only going to happen in the dense cen- 
tral region, where the baryons are going to form stars 
and subsequently remove the gas via their feedback. We 
simulate this effect by reducing mass of individual par- 
ticles in a central fraction of the system. The results in 
Figure [3] show that a more centralized mass loss leads 
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Fig. 4. — Evolution of the half-mass radius for intial conditions 
with a cusped profile {R^^ and R^^) or a turbulent velocity, re- 
spectively. 
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Fig. 5. — fcirci i-e. the potential in three selected cases. We plot 
the velocity in code units versus the radius in units of the final 
half-mass radius in the fiducial, undisturbed case. 

to a more extended asymptotic structure and makes the 
model more relevant for the formation of dwarf galaxies. 
As the dependance on the mass loss region is not sensi- 
tive, it merely demonstrates that an extended virialized 
system, as observed today, does not require mass loss to 
take place throughout the halo. 

Recent cosmological simulations (e.g. 

iDiemand fc KuhlenI I2008D show that dwarf galaxies 
form from the ins ide out, as in the clas sic secondary 
infall model (Fill more fc GoldreichI 119841 : iBertschingen 
|1985( ). To test the validity of our model in this context 
we investigate the collapse of a cusped profile. We 
set up two profiles, corresponding to p oc R^^ and 
p oc , respectively. As can be seen in Figure IH 

this leads to much smaller collapse factor. In fact, this 
is in better agreement w ith cosmological simulations 
(e.g. IDiemand et al.ll2007| ). as the simplicfication of a 
constant density is not applied there. However, the 
half-mass radius still is significantly puffed up in both 
cases with a cusped profile, similar to our simplified 
case. 

Another caveat in our initial conditions is the fact that 
the particles are at rest initially. To test this simplifica- 



tion we perform a run with a turbulent velocity. The 
velocity is set to correspond to half the circular velocity, 
i.e. the potential (see Equation [T] below) . Again, the 
maximum contraction is smaller than in the fiducial case 
(Figure HI . Still, in the case with mass loss, a significant 
change in the half-mass radius can be seen. We therefore 
conclude that the simplifications in our initial conditions 
do not affect the results strongly. 

In order to assess the change of the velocity dispersion, 
we evaluate the halo potential at a time t = 2tff, i.e. the 
circular velocity 



I'circ^y-^, (1) 

where G is the gravitational constant, and M is the mass 
inside the radius R. This approximation is adequate for 
small R where the system has already virialized. 

In Figure [SJ we plot Wcirc as a function of R/Rhm.s 
where i?hm,ff is the final half mass radius of the fiducial 
model (represented by the top black line) . We also com- 
pare the fiducial case and two models with a 10% mass 
loss where the mass loss is either confined to the cen- 
tral 20% (the bottom red line) or extended throughout 
the the entire halo (the middle blue line). These results 
show that mass loss significantly reduces the distribution 
of Wcirc at all radii to about 3/4 of the case without mass 
loss. 

Finally, the halos here all show a cored density profile , 
as has already been deduced by lAarseth et al.l ()1988D . 
This result is in agreement with the finding of most recent 
observational models whic h generally favor a cored rather 
than a cusped profile (e.g. lAmorisc o fc Evansii2012i) 

4. DISCUSSION 

The above models show, that a loss of 10% of the to- 
tal mass, from the central 20% of the total halo, at a 
time when the halo has contracted to 1/20 of its origi- 
nal size, is sufficient to reduce the velocity dispersion by 
^ 20 — 30% and to cause the halo to attain an extended 
asymptotic structure. 

In order to apply these simulated models in the con- 
text of Galactic dSphs, it is necessary to convert some 
of the dimensionless quantities into physical numbers. A 
straightforward quantity to do the mapping is the mass. 
By fitting the simulated and observed mass and velocity 
dispersion distribution, we can establish a characteristic 
mass, length scaling, and a complete normalization set. 

We apply this approach to fit the observed mass pro- 
file of the Fornax dSph. We adopt the model param- 
eters inferred from the ob served velocity dispersion by 
iJardel fc Gebhardtl (|2012l ). They determined its mass 
to be M300 = S-^toll X 10^ Mq inside R = 300 pc, 
M689 = 3.9t[! j^ X 10"^ Mq inside the projected half-light 
radius (i? = 689 pc) and M900 = 5.8to° x IO'^Mq inside 
the unprojected half-light radius {R = 900 pc). They 
also obtained an observed luminosity weighted line-of- 
sight velocity dispersion a = 11.3;ti'8 kms^^ 

We set the initial mass of Fornax dSph's entire halo to 
be 1.5 X 10^ Mq and choose a length scale which matches 
the observationally inferred mass distribution (Figure [6] 
left panel) . Then, we adopt their value of a and convert 

it to a circular velocity tidrc = VScr^ (|Wolf et al.ll2010l) 
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Fig . 6. — Mapping of the simulations to physical quantities 
ll20ll'l. Right: 

^circi i-^ - the potential, in three selected cases in 



IJardel fc GebhardtllSoH'l 

and plot that value at the position of the unprojectcd 
half-light radius (Figure |6] right panel). 

Although our choice of the initial total mass is some- 
what arbitrary, the combined constraint of the observed 
mass at various physical radii limits the degree of free- 
dom. There are several combinations of mass and radius 
normalization factors which can fit the data points in 
the left panel, but they all leave the match between the 
observations and our profiles in the right panel invari- 
ant. Therefore, we do the conversion with an reasonable 
choice of Mtot- 

Our choice of normalization also gives a timescale, tg = 
2.2 Gyr, which puts our simulations in a reasonable cos- 
mological context. The length scale is i?hm,fr = 9.2 kpc 
and the dark matter density is pdm ~ 3 x 10~^® g cm~'^. 
If we use a cosmological baryonic to dark matter ratio, 
the initial mass of baryonic mass would be ^ 2.55 x 
IO^Mq. Despite the loss of more than 50% of the initial 
baryonic mass (which corresponds to a loss of ~ 10% 
of the initial total mass), Fornax would still be able 
to have retained a sufficient amount of b aryonic mat- 
ter to account for its observed stellar mass (lMatealll998l : 
IJardel fc Gebhardt 20121) and some loss of baryonic mat- 
ter during the subsequent stellar and dynamical evolu- 
tion. 

Our best-fit model requires the mass loss to take place 
after the half mass radius has contracted by a factor of 
20. This collapse factor would reduce the instantaneous 
infall time scale to ~ 10 yr and the average density of the 
baryonic matter to ~ 10~^^ gcm~'^. These conditions are 
likely to lead to the onset of star bursts and consequently 
a strong feedback effect. A detailed investigation will be 
performed in future work. 

5. SUMMARY 

With the help of a simplified N-body model we are 
able to show a possible physical cause for the observed 
puff-up of satellite dSphs around the Milky Way. We 
propose that the internal star formation feedback during 
advanced stages of the first infall is sufficient to mod- 
ify the dark matter potential in order to explain their 



Left: Fitting the mass profile to observations by IJardel fc Gebhardti 
physical units in comparison with the observed value for Fornax (e.g. 

modified phase-space density. 

A natural next step is to use an SPH scheme to carry 
out a set of simulations which include a realistic treat- 
ment of the loss of baryonic matter. With the SEREN 
code, we can introduce a more realistic prescription for 
radiative or supernova driven loss of baryonic matter. 
For example, the impulsive mass loss prescription may 
be relaxed with a model for star formation which pro- 
ceeds on a time scale comparable to the dynamical time 
scale during the advanced stage of infall. Two types of 
particles may be used to represent gas and dark mat- 
ter. Star formation may mostly occur in the baryonic 
matter concentrated regions. We plan to include the 
Galactic potential which removes dark matter ventur- 
ing beyond the dSphs tidal radius. Recent simulations 
(e.g. lPawlik et al]|2012l) include some of these effects and 
show that the potential is indeed reduced in simulations 
with feedback (Brooks & Zolotov 2012). 

Here, we focus on the physical effect. During the 
first infall, the particles (dark matter and gas) are on 
near parabolic orbits. With the removal of a fraction 
of the mass, preferentially in the centre, the potential 
is changed sufficiently to make some particles orbits be- 
come unbound. We are able to show that e.g. for a loss of 
10% of the material in the central 20% of the (current) 
halo the observed lower value of the potential/the ve- 
locity dispersion can be explained within the framework 
of ACDM. The asymptotic dark matter distribution in 
these dSphs generally has a cored profile. 
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